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Context. Forsterite is one of the crystalline dust species that is often observed in protoplanetary disks and solar system comets. Being 



^0 I absent in the interstellar medium, it must be produced during the disk lifetime. It can therefore serve as a tracer of dust processing and 

disk evolution, which can lead to a better understanding of the physical processes occurring in the disk, and possibly planet fomation. 
However, the connection of these processes with the overall disk crystallinity remains unclear. 

Aims. We aim to characterize the forsterite abundance and spatial distribution in the disk of the Herbig Be star HD 100546, to 
Q I investigate if a connection exists with the large disk gap. 

5-H , Methods. We use a 2D radiative transfer code, MCMax, to model the circumstellar dust around HD 100546. We use VISIR Q-band 

ij^ • imaging to probe the outer disk geometry and mid-infrared features to model the spatial distribution of forsterite. The temperature- 

' dependent shape of the 69 /jm feature observed with Herschel** PACS is used as a critical tool to constrain this distribution. 

Results. We find a crystalline mass fraction of 40 - 60 %, located close to the disk wall between 13 and 20 AU, and possibly farther 
out at the disk surface. The forsterite is in thermal contact with the other dust species. We put an upper limit on the iron content of 
forsterite of 0.3 %. 

^ ' Conclusions. Optical depth effects play a key role in explaining the observed forsterite features, hiding waim forsterite from view at 

short wavelengths. The disk wall acts as a showcase: it displays a localized high abundance of forsterite, which gives rise to a high 
, observed crystallinity, while the overall mass fraction of forsterite is a factor of ten lower. 

oo , 

, Key words. Stars: individual: HD 100546 - Stars: pre-main sequence - planetaiy systems: protoplanetary disks - Radiative transfer - 

• . circumstellar matter 

o 

' 1. Introduction in the disk and settled to the mid-plane, depleting the disk of 

small, sub-micron sized particles. In addition, crystalline sili- 
Herbig Ae/Be stars are mtermediate-mass pre-mam-sequence ^^j^^ ^^^^^^^^ jl^^ mid-infrai-ed specti-a of many disks, 
. stars first described as a group by Herbig (1960). They are char- point i ng to substantial grain processing (e.g. Bouvvma^tlLl 
/\ ■ acterized by the presence of a circumstellar- gas and dust disk ^ Ivan Boekel et all '200.51. The recent discovery by direct 
^ , (e.g.| Waters&Waelkens| | 1998|) , which is the remnant of the star ima ging of e xo-planet s orbiting intermediate-mass stars (e.g. 
■ ■ ■ formation process. There is by now a wealth of observational Ma rois et al.l l2008: Kal as et alJl2008l:lLagrange et alJl20Tol) pro- 
evidence suggesting that these disks are the site of planet for- ^-^^^ ^^^^^^ ^^^^^^ jj^^ interpretation that disks surrounding 
matron. For instance, imaging at optical, mfrared and millime- ^^/g^ ^^^^^ ^j^^ .^^^^^ j^^^^ formation, 
ter wavele ngths shows the presence of disk gaps and/or inner 

holes (e.g. iGradv et al.l l200.5: Verhoelf et al. 201 1), pointing to However, many basic questions related to planet formation 
clearing of substantial parts of the disk. The composition of the "eed still to be clarified. For example: what is the main chan- 
dust in the disk strongly diff^ers from that of dust in the inter- "el for planet formation as a function of stellar mass? How does 
stellar medium. Large, millimeter sized dust grains have grown Pl^net formation aff^ect the composition of the gas and dust and 
its spatial distribution? Which observed properties of disks can 



Send offprint requests to: G.D.Mulders, e-mail: mulders@uva.nl be used as a signpost of planet formation, and can we derive an 
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** Herschel is an ESA space observatory with science instruments tems? In order to answer these questions, multi-wavelength ob- 

provided by European-led Principal Investigator consortia and with im- servations of the spatial structure and of the gas and dust compo- 

portant participation from NASA. sition of proto-planetary disks are needed. In this study, we focus 
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on the spatial distribution of crystalline silicates in the Herbig 
Be star HD 100546. This star has been observed with Herschel 
PACS as part of the DIGIT Q open time key pr ogramme that aims 
to address these questions (iSturm et alJbOlOl) . 

Crystalline silic ates are not found in interstellar space 
(iKemper et alJ l2004l) but are abundant in a significant fraction 
of pr oto-planetary disks (e.g. JOlofsson et 311120091: [juhasz et alJ 
12010 ). They must therefore have been formed in situ, and trace 
the thermal and chemical history of the grains in the disk. 
They can form by thermal annealing of amorphous silicates 
(iFabian et al.l l2000l) . or by direct condensation from the gas 
phase. Both processes require high temperatures, respectively 
above the glass temperature (~1000 K) or near the dust evapo- 
ration temperature (~1500 K). These temperatures are typically 
found close to the star. 

Crystals are also observed much farther away from 
the star than expected based on their temperature (e.g., 
lOlofsson et al.l 12010) and are abundant in solar system 
comets, requiring eflicient rad ial mixing o r local production i n 
shocks (i Harker & D esch'2002'). colHsions jMorlok et al.ll2010l) . 
parent-body processing (Huss et al. 200 Ij) or stellar outbursts 
(I Abraham et alJ l2009h . Additionally, crysta ls can also be de- 
stroyed by stellar winds (iGlauser et al.ll2009b . However, the cor- 
relation of crystallinity with other stellar parameters such as age 
and luminosity in a large sample remains unclear (e.g. Watson 
[2009; Oliveira et al. 2010; Juhasz et al. 2010) . Detailed studies 
of individual disks are therefore needed in which spatially re- 
solved information and spectroscopy are combined in order to 
establish the spatial distribution, abundance and chemical com- 
position of crystalline silicates in relation to the overall disk ge- 
ometry. 

In this paper we focus on HD 100546, a Herbig Be star of 
spectral type B9.5Vne with one of the strongest observed crys- 
talline dust features (Malfait et al. 1998). Apart from the un- 
usual mineralogy, the Spectral Energy Distribution (SED) of HD 
100546 is characterized by a large mid-infrared excess and a 
modest near-infrared excess. This SED is explained by a disk 
gap, which creates a frontally illuminated wall at the far edge of 
the gap with a temperature of ~ 200 K (Bouwman et al. 2003). 
The existence of the gap has been confirmed by spatially re- 
volved s pectroscopy of the gas as well, placing the wall at 13 
AU re.g. lGradv et all 120051) . The small near-infrared excess is 
caused by a tenuous inner disk (Benistv et al. 2010). 

The m id-infrared spectrum of HD 10 0546 has been observed 
with ISO (iMalfait et all 19981). T IMMI2 llvan Boekel et alJ2005h 
and Spitzer (Juhasz et al."2010|). The ISO observations revealed 
very pronounced crystalli ne silicate feat ures, similar to those of 
comet Hale-Bopp (iCrovisier et al.l[T997l) . Comparison to labora- 
tory spectra indicated that the features between 1 1 and 30 fim 
can be de scribed by a dust p opulation with a single temperature 
of 210 K (IMalfait etal J 1998 ). coinciding with the disk wall. The 
longer wavelength features (33 and 69 /vm) require an additional 
cold component (~50-70 K) to explain the feature strength. 

The ISO spectru m was analysed in more depth by 
iBouwman et alj (l2003h . who study the forsterite abundance and 
spatial distribution with a radiative transfer model that assumes 
the dust is optically thin. They find a crystalline mass fraction 
that increases from 2% in the inner disk to 19% in the outer disk 
- which is inconsistent with radial mixing models that predict a 
decreasing fraction with radius. 

The TIMMI2 spectrum was analysed using a similar opti- 
cally thin model, but without taking into account the geometry, 

' Dust, Ice and Gas In Time 



bv Ivan Boekel et aT] (l2005h . who find a forsterite mass fraction 
of 5.0% in the 10 ;um range. The forsterite features have been ob- 
served more recently with Spitzer by Juhasz et al. (2010). They 
find a forsterite mass fraction of 5.6% from the 10 fim spectral 
range and 5.1% from the 20 fim spectral range. 

In addition to the strength, the shape and central wavelength 
of the features also contain information on the temperatur e of the 
emitting dust, as these are temperature-dependent (Koi ke et al.l 
2006). This is especially true for the 69 fim feature where the 
feature broad ens and sh ifts to longer wavelengths at higher tem- 
peratures (Bo wev etai] l2002; Suto et al. 2006, Fig.lCTTl. 

The 69 fim feature of HD 100546 has been observed in great 
detail with Herschel PACS a s part of the DIGIT open time key 
programme, and analysed by ISturm et aP (1201 Oh using an opti- 
cally thin approach. The observed central wavelength of the fea- 
ture is 69.2 fim, whereas cold (50 K ), iron-free forsteriteQ has a 
peak at 69.0 fim. ISturm et al.l (1201 Ol) propose two different sce- 
narios for the required wavelength shift of 0.2 fim: 

In the first scenario, warm iron-free forsterite creates the ob- 
served shape. Using a weighted sum of the opacities, they find 
that pure forsterite with a temperature of 200- 1 50 K has a peak at 
69.2 fim and matches the feature shape. This scenario is consis- 
tent with the inferred temperature of the disk wall, but strongly 
overpredicts the features at shorter wavelengths. It is therefore 
inconsistent with the ISO and Spitzer data, unless optical depth 
effects play a significant role. 

In the second scenario. lSturm et al1(l2010h proposed an alter- 
native explanation: colder forsterite, with a temperature around 
50-70 K creates the 69 fim feature, and an admixture of a few 
percent iron shifts the peak longwards to 69.2 fim. Such a model 
does not overpredict the sh ort wavelength featu res, consistent 
with the earliest model from IMalfait et al.l (Il998l) . Optical con- 
stants for forsterite with a few percent iron have not (yet) been 
measured in the lab, and a compositional fit could not be made. 
The inferred wavelength shift follows from in terpolation of mea - 
surements with 0% and 10% iron content (see lSturm et al.l2010l) . 

The two scenarios lead not only to a diff'erent chemical com- 
position of the forsterite, but also to a different location of emis- 
sion within the disk. Iron-free forsterite originates from the disk 
wall, whereas iron-containing forsterite comes from colder ma- 
terial farther out - possibly near the disk midplane. Because 
forsterite has a large number of features over a broad wavelength 
range, regions with different optical depths are probed at dif- 
ferent wavelengths, providing extra diagnostic power To shed 
more light on the chemical composition and spatial distribution 
of forsterite in the disk, we need to explore how feature strengths 
and shapes over the entire wavelength range are influenced by 
these optical depth efifects. 

To do this, we need to use a 2D radiative transfer code, de- 
scribed in Section |3] to construct a physical and mineralogical 
model of the dust surrounding HD 100546. We construct a struc- 
tural hydrostatic model in section |4] which is consistent with a 
new VISIR mid-infrared image (appendix lAli. In section |5] we 
use this structural model to study the spatial distribution of crys- 
talline forsterite, by fitting the forsterite feature strengths over a 

- Ciystalline silicates with olivine stoichiometry are called forsterite 
in the case they contain only magnesium (Mg2Si04), and fayalite if 
they contain only iron (Fe2Si04). In the case they contain both iron 
and magnesium (Fe,.Mg2_.vSi04), they are called olivine. Because of 
the extremely low iron content of the crystalline olivine (x < 0.02) used 
in this paper, and to avoid confusion with amorphous olivines, we will 
refer to them as forsterite as well. We will use the denomination iron- 
free / iron-containing when it is necessary to discriminate between both 
species. 
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broad wavelength range and in particular the shape of the 69 fim 
feature. 

We will discuss our results in the context of previous work 
in section |6l as well as address the question of the origin of the 
forsterite and correlation with the disk gap. We will summarize 
our conclusions in section |7] 



2. Observations 

2.1. SED and spectroscopy 

To construct the observed SED of HD 1 00546 (see Figu re |2]i, 
we use photometric data from van den An cker et al.l(ll997b com- 
bined with the ISO spectrum fromMalfait et al.l (Il998l). We do 
not u se the Spitzer (lJuhasz et al]|2010h and PACS lSturm et al.l 
12010 ) spectra for constructing the SED because the ISO spectra 
encompasses the wavelength range of both, and the current flux 
calibration of the PACS spectrum is not yet better than that of 
ISO. There is no significant offset between the ISO and Spitzer 
spectra that could influence our analysis. We do make use of the 
increased sensitivity of PACS and Spitzer for calculating the in- 
tegrated strength of the forsterite features. 



Dust mineralogy 



3. Radiative transfer and dust model 

3.1. Radiative transfer model 

The dust radiativ e transfer code used in this paper is MCMax 
(iJVIin et alJl2009h . a 2D Monte Carlo code. It is based on the im - 
mediate re-emission procedure from Biorkman & WoodI (12001 1, 
combined with the method of continuous absorption by Lucvl 
(119991) . The code has been benchmarked against oth er radiative 
transf er codes for modelling protoplanetary disks (iPinte et al.l 
l2009l) and has been successfully applied f or modelling spatiall y 
and spectrally resolved observations (e.g. Verhoeff et al 



I2010h . 



The radial grid around the inner radius and disk wall is refined to 
sample the optical depth logarithmically. The SED, images and 
flux contributions (see section |53] | are calculated by integrating 
the formal solution to the equation of radiative transfer by ray- 
tracing. 

In addition to radiative transfer, the code explicitly solves for 
the hydrostatic vertical structure of the disk, with the implicit as- 
sumption that the gas temperature is set by the dust temperatureQ. 
Therefore it only needs the radial distribution of dust to construct 
the structural model. For the radiative transfer that sets the dust 
temperature - and therefore its vertical structure - also the stellar 
and dust properties are required. 

The radial distribution of the dust in the disk is defined by 
the following parameters: the inner and outer radius (Ri,, and 
Rout) as well as the gap location Rgap.in and Rgap.out; the total 
dust mass (Mdust) and its distribution across the disk, the surface 
density profile (SDP); and finally an additional dust depletion 
factor for the inner disk finnei- The stellar spectrum is described 
by a Kurucz model with effective temperature Teff, luminosity 
mass, and is set at a distance d. These parameters are 
summarized in Table |2] 

Although th e stella r parameters have recently been updated 
by iTatulli et"al] (1201 lb . our fit parameters provide an equally 
good fit to the stellar photosphere. Adopting the lower luminos- 
ity of 26 Lq would have only a small impact on the model, as it 



Dust species 




reference 




13.8% MgFeSi04 


[1] 


Amorphous 


38.3% Mg2Si04 


[2] 


silicates 


42.9% MgSiO, 


[1] 




1.8% NaAlSiaOe 


[3] 


Carbon 


C 


[4] 


Forsterite 


Mg2Si04 


[5] 



Table 1. Dust composition used in this paper Opacities are 
calculated from the optical constants using a size distribution 
from 0.1 to 1.5 jum proportional to a"^^. The shape of the par- 
ticles is DHS, with fmax= 0.7 except for crystalline forsterite, 
which has fmax= l O- Refe rences to the op tical con stants : [1] 
Dorschn eret"an (11995): [2] Henning & Stognienko (199 6|); [3] 
■MutschkeetalJ([l998.) : r41 Preibisch et al..(.1993,) : r51 .Suto et alJ 
120061) 



requires a higher continuum opacity in the optical to achieve the 
same dust temperatures and fit the SED. 

3.2. Dust model 

We compute the optical properti es of our grains us ing a distri- 
bution of hollow spheres (DHS, iMin et alJ l2005ah . which has 
been shown to rep resent mid-infrared feature shapes very well 
(lJuhasz et al.ll2010l) . 

We will restrict ourselves to a model with three dust com- 
ponents that dominate the opacity from the optical to the far- 
infrared: amorphous silicates, crystalline forsterite and carbon. 
Amorphous silicates produce the 10 and 20 fim emission fea- 
tures which are very prominent in HD 100546. We will use 
the composition and optical constants of i nterstellar s ilicate s 
towards the galactic center as derived by iMin et al" ] dloOTl) . 
which are summarized in table [1] Crystalline forsterite dom- 
inates the mid infrared opacity only at wavelengths where it 
has strong features, most notably at 11.3, 23.5, 33.5 and 6 9 
yum. We will use the optical constants from ISuto et all (l2006l) . 
as these have been measured at different temperatures. At opti- 
cal and near-infrared wavelengths, the opacity is dominated by 
continuum sources without spectral features in the mid and far- 
infrared. For this component we use carbon, with optical con- 
stants from Pr eibisch et al.l (Il993b . Note that probably not all of 
the continuum-opacity sources are amorphous carbon, and that 
also different materials with similar optical properties could con- 
tribute, for example metallic iron. 

Dust grains in protoplanetary disks are most likely in the 
form of mixed aggregates. Lacking a good effective medium the- 
ory for DHS, we calculate the opacities for every dust species 
from the optical constants as if they were separate particles, and 
then force them to be in thermal contact. Although this approach 
does not have as good a physical basis as an effective medium 
theory, it does allow an accurate representation of feature shapes. 
In addition, opacities do not have to be recalculated when the 
composition changes, making it easier to conduct large parame- 
ter studies and make models with spatially varying dust compo- 
sitions, including temperature dependent opacities. The implica- 
tions of the dust grains being in thermal contact will be discussed 
in sectionl6^ 



The gas scale height is therefore not a free parameter in our model. 



The value derived bv lvan den Ancker et al.llll997h is L,=33 Lq, but 
L,=36 Lq, which is within the Icr error, better fits the available photom- 
etry. 
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Model parameters 



Parameter 


Value 


reference 


Agff 1^1 


10500 


n 1 


T n 1 


36+ 


n 1 




2.4 


ni 

J 


d fncl 


103 


ni 

J 


R TATTl 


0.25 




Rexn fAUl 
exp L ^ J 


350 


[3] 


Rou( [AUl 

^^OIU L ^ J 


1000 


[4] 


R„,ni„ [AU] 

gap, in L ^ J 


4 


[2] 


-Rgap,out [ A^--^] 


13 


[51 

J 


Mhm,, TMnl 


1 X lO-* + 






r-' + 




f. 

Ainner 


200+ 




amin[/Jm ] 


0.1 


[6] 


amax[j"m] 


1.5 


[6] 


carbon fraction [%] 


5 + 




Shape 


in-egular (DHS) 


[6] 


in 


42 


[7] 


PA 


145 


[7] 



Table 2. Model parameters. All parameters with a dagger (') 
are (re)fitted. Refere nces: [1] van den Ancker et al. (1997); [2^ 
iBenistv et alJ (l2010l): [3 1 jPanic et al. (2010); [4] Ardilaet; 
(I2007h: [51 iGrady et al.l EoOS.) : [6] Juhasz et al. (.20101) 7" 
iPantin et all hoOOl) 



[2] 
[7] 



Other dust properties - grain size range and composition - 
have to be constrained from the observations and are discussed 
in section |431 



4. Disk structure 

Before we can model the spatial distribution of forsterite, we 
need to construct a structural model for the disk that fits the avail- 
able observations. The disk of HD 100546 has been well-studied, 
and the disk geometry is already well-constrained. In the next 
section we will describe these observational constraints. In sec- 
tion l4.3l we fit the SED to fill in most of the remaining parameters 
of our disk model. The only parameter important for this study 
that can not be well constrained from the SED is the surface 
density profile (SDP). We will assume the same SDP as previ- 
ous modeling attempts, which is consistent with VISIR imaging 
(see appendix IA.3I ). Fig [T] shows a Q band image of the final 
disk model, illustrating its geometry with an outer disk (Section 
14.1.1b . gap (Section l4T2t and inner disk (Section l4T3l l. 

4.1. Observational constraints 
4.1.1. Outer disk 

The (out er) disk of HD 10054 6 was first resolved in scattered 
light by Aug ereau et al.l (|2001[). and has since then been image d 
at different wavelengths (e.g. lLiu et al.l2003LlArdila et al.ll2007h . 
Its position angle and inclination have been measured using dif- 
ferent techniques, indicating an inclination between 40 and 50 
degrees from face-on, and a position angle east of north between 
130 and 145 degrees (e.g. lPantin et al.ll2000l: lArdila et al.l2007h . 
We adopt the inclination of 42 degrees and a position angle of 
145 degrees. 

The outer radius of the dust disk is difficult to determine. 
Scattered light is observed up to at least ten arcseconds (1000 
AU at a distance of 103 pc. lArdila et al.ll2007l) . yet it is not clear 
whether this belongs to the disk or to a remnant envelope, and the 




Fig. 1. 18.7 fim model image using our final model of the inner 
regions of HD 100546. Note that the inner disk casts a clear 
shadow on the outer disk wall. The field of view of the image 
is 0.5" by 0.5", or 51.5 AU at a distance of 103 pc. North is up 
and East is left. 



o 




Fig. 2. SED of HD 100546 (diamonds). The grey fine is the ISO 
spectrum. Overplotted is the best-fit disk model with forsterite 
(R20, solid line) and without forsterite (dotted line), both with- 
out PAHs. 



dust disk could be significantly smaller (e.g. lArdila et al.ll2007h . 
Modelling of rotational CO lines are co nsistent with a ga s-rich 
disk with an outer radius of 400 AU (IPanic et al.' '2010'). We 
therefore adopt an outer radius of 1000 AU, but with an expo- 
nential cutoff that sets in at Rexp= 350 AU (iHughes et al.ll2008h . 

The total dust mass of the outer d i sk is 7 .2 X 10""* M©, as 
has been determined by Henning et al.l (Il998l) from 1 . 3 mm mea- 
surements, which probe mostly big grains in the disk midplane. 
The mass distribution in the outer disk is not well constrained. 
Vari ous authors have adopted an SDP of £ (r) - r^P with p=l 
(e.g. iPanic et al.ll2010l iBenistv et alJl2010l) . a value commonly 
used for protoplanetary disks. We will also adopt this value, 
which is consistent with our VISIR image in appendix lA.3l 

4.1.2. Gap size 

The SED of HD 100546 (Fig. ^ is dominated by a relatively 
small near-infrared excess between 2 and 8 fim, followed by a 
steep rise towards the mid-infrared at 10 //m. This SED has been 
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attributed to a disk gap around 10 AU - possibly cleared out by 
a proto-planet - which creates a frontally illuminated wall at the 
far edge of the gap (Bouwman et al. 2003). This wall has a tem- 
perature of 200 K, and is responsible for the huge mid-infrared 
excess compared to other Herbig Ae/Be stars. The gap mani- 
fests itself in mid-infrared imaging, where the source is much 
larger than other Herbig Ae/Be stars, and its size is independent 
of wavelength ( Liu et al, 2003) . 

The gap is also seen in the distribution of the circum- 
stellar g as (e.g. Acke & van de n Anckerl l2006i iBrittain et al.l 
[2009 v an der Plas et al. 2009), and tight constraints have been 
put on the size of the disk ga p using spatially resolved spec- 
troscopy bv^ Gradv et al.l (l2005b . We adopt their value for the gap 



outer radius of R, 
wall. 



gap,out" 



= 13 AU, which is the location of the disk 



4.1 .3. Inner disk 

The inner disk of HD 100546 causes a near-infrared excess 
that is relatively small compared to other Herbig Ae/Be stars 
(Bouwman et al. 2003), indicati ng that the dust is depleted 
(iDominik et al., ,2003) or settled. iBenistv et aP (1201 >) have re- 
solved the disk using near-infrared interferometry, constraining 
the inner radius to Rin=0.25 AU. The authors also find that the 
inner disk is very tenuous, and has a low dust mass of only 
3 X 10^'" Mq. The AMBER observations cannot resolve the 
outer radius, but they put it at Rgap,in=4 AU. We adopt these val- 
ues for the inner and outer radius. 



4.2. Mineralogy 

Before determining the remaining model parameters using SED 
fitting we need to specify our dust composition. Although the pe- 
culiar forsterite mineralogy of HD 100546 is the main focus of 
this paper, forsterite itself plays only a minor role in determining 
the disk structure and fitting the SED. The main reason for this 
is that it has a very low opacity in the optical compared to other 
dust species such as amorphous silicates and carbon. The latter 
species therefore determine the heating/cooling balance of the 
dust, and set the disk's temperature and vertical density struc- 
ture. Leaving out the forsterite does not affect the shape of the 
modeled Q band radial profile, and changes to the overall SED 
are small (Figure|2}. The shape of the 10 micron feature is better 
fitted with forsterite included, but this does not affect the overall 
SED shape (see Figure |2]l. We will therefore constrain the disk 
geometry in this section using only the other two components, 
amorphous silicates and carbon. By doing so, the derived ge- 
ometry is independent of the forsterite mineralogy, and we can 
study the abundance and spatial distribution of forsterite in sec- 
tion |5] without refitting the SED. 

The dust composition is determined by the mass fraction 
of carbon, and f ollows fro m our SED fi t in se ction 14.31 Both 
iBouwman et alJ (l2003h and lJuhasz et all (l2010l) show that the 
small grains responsible for the mid-infrared dust emission fea- 
tures in HD 100546 cannot be larger than a few //m, though 
larger grains could be present at higher optical depths not probed 
by Spitzer. The small grains are most important for setting the 
disk dens ity and temperatur e structure because of their high 
opacities (iMeiier et al.ll2008 \ and we adopt a grain size distri- 
bution between 0.1 and 1.5 yum for all grains. We do not fit the 
SED in the (sub)mm because it probes a separate population of 
large, millimeter-sized grains (iBenistv et al.l2010l) . which do not 



influence the strength of mid-infrared features. A summary of 
the mineralogy can be found in table [T] 



4.3. SED fitting 

With the disks radii (Rin, Rgap.in, Rgap.out, Rout) and opacities con- 
strained, the next step is to constrain the dust mass (Mdisk), inner 
disk depletion (finnei) and carbon fraction by SED fitting. The 
SED is relatively insensitive to the SDP, but some limits can be 
placed on it using spatially resolved VISIR imaging (See ap- 
pendix |A]l. For now, we will adopt p - I. The total dust mass 
for HD 100546 has been measured at millimeter wavelengths 
(iHenning et al.|[T998h . However this mass probes mostly the big, 
milUmeter sized grains in the midplane, while we are interested 
in the dust mass in small micron-sized grains responsible for 
mid-infrared and far-infrared features. We therefore have to re- 
fit the dust mass to the optically thin part of the SED at ~ 
60-200 //m. We find a dust mass of 10~^ Mp, con s istent with 
the dust mas s in s mall grains from iDominik et al.l (l2003l) and 
IBenistv et all (^0). 

We proceed by fitting the inner disk (0.25-4 AU), which has 
been done before w ith a similar 2D radiative transfer code by 
IBenistv et al.l (12010 ^. They use a parametrized vertical structure 
with a fixed scale height and flaring exponent. In our model, the 
vertical structure is computed under the assumption of vertical 
hydrostatic equilibrium. Therefore we can not freely adjust these 
parameters, so instead we fit the inner disk SED by varying its 
dust mass. We do this by assuming the mass distribution in the 
inner disk follows the same power law as the outer disk {p - \) 
- as would be the case for a disk without a gap - and reduce 
the surface densit}0 to lower the opti cal depth and fit the near- 
infrared flux in the same way as Ben istv et al.l (l2010h . We fit a 
dust mas s of 6 x 10"^ M©, consistent with the value found 
by Benist" vet al.l (|2010)- This corresponds to a surface density 
normalization of Siau = 2.3 x 10~^ g/cm^ . This is a reduction 
of a factor finnei— 200 with respect to an extrapolation inward of 
the SDP of the outer disk, which has Eiau = 0.45 g/cm^ . 

The modelled inner disk is not completely optically thin 
in the radial direction at optical wavelengths, as can be in- 
ferred from the shadow seen on the disk wall in figure [1] Even 
when we take into a ccount that the star is not a point source 
(EspainatltiD|2010!), the shadow reaches up to a height above 
the midplane of around ~ 1 AU at the disk wall. The shadow 
only disappears when the inner disk has an incre ased scale height 
such that it effectiv ely becomes a halo (Mu lders et al.l I2010I: 
'Verhoeffet al."201l'), as has been suggested for HD 100546 by 
Vinkovic et al. (2006). 

Even though the wall at 13 AU is partly shadowed, the il- 
luminated part still intercepts a significant fraction of the stel- 
lar light. It is frontally illuminated by the star, and has a tem- 
perature of 200 K. Such a structure fits the mid-infrared part 
of the SED without any increase in scale height beyond hydro- 
static equilibrium - which is often needed to a ccount for extrem e 
near-infrared excesses in Herbig Ae/Be stars ("Acke et ani2009l) . 
The steep rise in the SED between 8 and 10 yum is caused by 
this wall, but it also contains information on the dust composi- 
tion. At 8 //m the opacity is dominated by carbon, at 10 jum by 
amorphous silicates. From the 10/8 //m flux ratio, we estimate a 
carbon/silicate mass ratio of 5%. Increasing the carbon fraction 
'fills' the gap in the SED at 8 //m: with 10% carbon the flux at 



' Reducing the surface density has a similar effect on the near- 
infrared flux as reducing the scale height or increasing the flaring ex- 
ponent. 
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Fig. 3. The temperature structure of the outer disk of HD 1 00546. 
Overplotted are the radial t-1 surface in the optical (solid thick 
line) as well the vertical t=1 surfaces of the continuum at 24, 
33 and 69 jum (dotted, dashed, dot-dashed). Note that the verti- 
cal T-l surface is calculated from above, meaning that regions 
with higher optical depth lie below it. The 69 fim surface lies be- 
low the disk midplane, meaning that the disk is close to getting 
optically thin at these wavelengths. 
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Fig. 4. Continuum subtracted spectra comparing the optically 
thick (solid black) and optically thin (light grey area) approach 
using only warm forsterite (model RW30). Models are fitted to 
the 69 jum feature shape. Overplotted are the Spitzer spectrum 
between 7 and 35 fj.m and the PACS spectrum between 65 and 
73 jum (dark grey lines). A more detailed version comparing the 
optically thick model and observations can also be found in fig- 
ure |5] panel RW30. 



this wavelength is already significantly ov erpredicted. Because 
carbon combusts (see, e.g., Lee et al. 2010) at the same temper- 
atures needed to anneal silicates ( Gail 200 1). the relatively low 
carbon abundance could be linked to the high crystallinity. 

4.4. Final model 

Model parameters of the final disk model used for the anal- 
ysis of the forsterite features can be found in table |2] Note 
that this model is in hydrostatic equilibrium, and therefore does 
not have a predefined scale height. The calculated scale height 
increases with height above the midplane due to the temper- 
ature gradient, but can be fitted with a power law in radius 
(H= HiooAU * (r/lOOAU)'') at both the surface and midplane for 
comparison. The outer disk is well described by a flaring expo- 
nent of yS = 1 .4 and a scale height that increases from H iooau =13 
AU in the midplane to Hiooau=21 AU at the disk surface and in 
the disk wall. The inner disk is fitted withjS =1.3 and Hiooau=13 
AU in the midplane to Hiooau=25 AU in the disk rim. The mid- 
plane temperature follows T = ~150 * (r/AU)""'^^ for the outer 
disk, and is a factor of 2 higher in the inner disk. 

5. Disk mineralogical model: Forsterite 

In this section we will study the abundance and spatial distribu- 
tion of forsterite in the disk of HD 100546, using the disk struc- 
tural model described in the previous section. The aim is to dis- 
tinguish betw een the two s c enario s for the forsterite composition 
described by ISturm etaH (1201 Ob (warm iron-free forsterite or 
cold iron-containing forsterite), which have also been discussed 
in the introduction. 

5.1. Optical depth effects 

We explore in detail the effects of optical depth and the gapped 
disk geometry of HD 100546 on the observed feature strengths 
of forsterite. Throughout the paper, we will assume that the 



forsterite is in thermal contact with the other dust species (see 
section |4~2] ). Fig [3] shows the temperature structure in the outer 
disk from the model derived in section l4~4l Fig|3]also shows the 
vertical t=1 surfaces - as seen from the observer to the emission 
point - of the continuum at 24, 33 and 69 fim. These surfaces 
mark the boundaries below which the dust can no longer con- 
tribute substantially to the feature at that particular wavelength. 

At 69 fim, the vertical t=1 surface lies below the disk mid- 
plane (i.e., the disk becomes optically thin), allowing most of the 
disk to contribute to the 69 fim feature. At shorter wavelengths, 
these surfaces lie above the disk midplane out to radii of 50-100 
AU. For the 150 and 200 K components, a significant part lies 
below the t= 1 surface, especially when taking into account that 
the density increases towards the midplane. Most of the warm 
forsterite that can contribute to the 69 /im feature is therefore 
hidden from view at shorter wavelengths. 

To illustrate the influence of these optical depth effects on 
the measured feature strengths of forsterite, we compare the 
continuum subtracted spectra (see appendix |B]i of our model 
(described in secti on I5.2.4l i with the optically thin model from 
ISturm etaTI (1201 Oh in Figure|4] Both models are fitted to the fea- 
ture shape and strength at 69 fim, and contain only 150 K and 
200 K forsterite. The optically thin approach clearly overpre- 
dicts the strength of the short wavelength features in the 10 to 40 
fim range, while the 2D radiative transfer does not overpredict 
these because the bulk of the warm forsterite is hidden at such 
short wavelengths. 

The gapped geometry of HD 100546 plays a key role here: 
because part of the wall at 13 AU is directly illuminated by the 
central star, it heats regions close to the midplane to temperatures 
of 150 to 200 K. In a disk without a gap, these regions would be 
much colder and these high temperature regions would only be 
found in the midplane at much smaller radii (~4 AU). There, the 
surface density would be high enough also for the t-1 surface 
at 69 /im to be above the midplane. Without a gap, the optical 
depth effects would not be so pronounced as the hot forsterite 
would be hidden from view at 69 //m as well. 
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This first analysis shows that while optically thin models 
with iron-free forsterite overpredict feature strengths at shorter 
wavelengths, a 2D radiative transfer model that takes optical 
depth effects into account does not. Additionally, none of the 
longer wavelength features (33, 69 fim) require an additional 
cold (<150 K) component to fit the feature strengths. This indi- 
cates that the shape of the 69 fim feature can be fitted with only 
warm forsteri te, suggestin g an iron-free composition (the first 
scenario from lSturm et al.l uOlO). We will therefore first pursue 
a disk model that explains all feature strengths and shapes with 
iron-free forsterite. Afterwards, we will discuss models that in- 
clude a fraction of iron. We describe how we implement the 
temperature-dependent opacities of iron-free forsterite into our 
radiative transfer code in Appendix ICl 

5.2. Forsterite spatial distribution 

To determine the spatial distribution of the forsterite in the outer 
disk, we use the disk structural model from section l4~4l We re- 
place a small fraction of the amorphous silicates by forsterite, 
and compare continuum subtracted model spectra with the data. 
All discussed models can be found in Fig. |5] These models are 
fitted to the integrated flux at 69 fim, which is ~9 Jy fim. 

We will show that models with constant abundance through- 
out the outer disk severely underpredict short wavelength fea- 
tures due to optical depth effects (model GO). We therefore 
run an extensive parameter study where the radial distribu- 
tion of forsterite peaks towards the disk wall. We consider two 
parametrizations: an abundance gradient (model series G), pre- 
dicted by radial mixing models (Gail 2001), and a (narrow) ring 
of constant abundance (model series R), which makes it easier to 
control the minimum temperature of the forsterite. In addition, 
we also ran a series of models where the forsterite is only located 
in the disk wall and at the disk surface (model RW), and not in 
the midplane. Below we discuss the different models. 

5.2.1 . Constant abundance throughout the disk 

In contrast with the result from iBouwman et al. I (l200l . models 
with a constant forsterite abundance throughout the outer disk 
consistently fail to reproduce all observed feature strengths. This 
failure is due to optical depth effects, which suppress short wave- 
length features (Section FS.lb . A very low abundance of 2 % fits 
the integrated 69 fim flux (Fig.|5] model GO), but underpredicts 
shorter wavelength features, most notably shortwards of 33 fim. 
An abundance increase of at least a factor of 10 is necessary in 
the disk wall to fit the flux at 1 1 fim. To achieve this without 
overpredicting the 69 fim feature by the same factor, we need to 
concentrate the forsterite towards the disk wall. 



5.2.2. Abundance gradient 

First, we concentrate the forsterite towards the disk wall by in- 
troducing an abundance gradient that gradually decreases out- 
ward. Radial mixing models predict that if radial mixing trans- 
ports crystalline material fr om the inne r to the outer disk, it fol- 
lows a smooth distribution (iGaifcOOlh . which is well described 
by a power law (IStevensonlll990h . We parametrize the forsterite 
abundance as /(r) - /hau ir[AU] / 13 AUy^ , where free fitting 
parameters are the abundance in the disk wall (f^Au) and the 
steepness of the gradient (^). The grid scanned for these parame- 
ters is summarized in table|3] We start by assuming the forsterite 



Forsterite spatial distribution 
Parameter Best-fit Range explored 

abundance [%] 40 {2, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90) 

rn„B [AU] 20 (14, 15, 16, 17, 18, 19, 20, 25, 

30, 35, 40, 60, 80, 100, 400, 1000) 
_f 2 {0, 0.5, 1.0, 2.0,3.0, 4.0) 

Table 3. Forsterite abundances in the disk. Best-fit and ranges 
explored. Abundance refers to abundance at 13 AU; outwards 
the abundance is either constant up to rring, beyond which it is 
zero, or it drops off gradually as ( '^^au 



is well-mixed in the vertical direction, i.e. there is no vertical 
abundance gradient. 

With a gradient of ^=1, we can increase the forsterite abun- 
dance in the wall to 10% without overpredicting the 69 fim flux. 
This is still not high enough to explain the short wavelength fea- 
tures, which originate mainly in the disk wall (Gl). Steeper gra- 
dients of ^=2 or 3 concentrate the forsterite closer to the disk 
wall, and can reach higher abundances for the same 69 fim flux. 
They reach abundances of respectively 30% (G2) and 40% (G3), 
which is sufficient to fit the feature strengths at shorter wave- 
lengths. With even steeper gradients, the abundance drops off so 
rapidly even in the wall that short wavelength features are again 
underpredicted (not plotted). 

Although models G2 and G3 can explain all feature 
strengths, the shape of the 69 fim feature does not fit the observa- 
tions well. The peak is located at 69.0 fim, indicating that cold, 
100 K forsterite dominates the peak location. Although an in- 
creasing contribution at 69.2 fim is seen with steeper gradients, it 
never dominates the feature shape. These solutions therefore re- 
quire some iron to shift the feature towards longer wavelengths. 
However, even then the feature shapes do not match very well 
(Figure |5]l. 

The reason why the cold forsterite is so prominent at 69 fim 
is that even with a steep gradient a few percent forsterite remains 
far out in the disk. Such a low abundance provides a significant 
contribution at 69 fim, while it does not add a lot of flux to the 
short wavelength features (see also model GO). We will therefore 
focus on a series of models with a forsterite abundance that is 
non-zero only near the disk wall. 

5.2.3. Abundance ring 

To prevent a low percentage of cold forsterite in the outer 
disk from dominating the shape at 69 fim, we choose a second 
parametrization where the forsterite is located in a narrow ring 
starting at the disk wall. The abundance of forsterite is constant 
within the ring, between 13 and rrfng AU, and is zero outside this 
ring. Both the abundance in the ring and its outer radius r^ng are 
free fitting parameters. The explored parameter space is summa- 
rized in Table [3] a selection of models is displayed in Figure |5] 

Because the midplane is close to optically thin at 69 fim, 
avoiding a strong contribution from cold (50-100 K) forsterite 
requires a narrow ring of only a few AU. A ring with outer ra- 
dius of 14 AU (R14) requires a high abundance of 60% to fit the 
integrated feature strength at 69 fim. At such high abundances, 
features become optically thick and flatten off, producing a 69 
fim feature that is much too broad. Additionally, the short wave- 
length features flatten off, and become too weak compared to the 
observations. 

Lower abundances require broader rings to produce enough 
emission at all wavelengths. Rings with abundances of 40 and 
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Fig. 5. Continuum subtracted spectra comparing the models (solid black) and Spitzer and PACS observations (grey). Presented 

models have integrated feature strength between 67.5 and 71.0 //m of 8.5 ± 1.5 Jy. Offsets between models are 50 and 15 Jy for the 

left and right panel respectively. The cartoons on the right show the location of the forsterite in the outer disk, not to scale. A higher 

intensity represents a higher abundance. Roman symbols present models with: 

GO - 2% forsterite throughout the entire disk (gradient with r"*^); 

Gl - gradient with 10% forsterite at 13 AU and r"^; 

G2 - gradient with 30% forsterite at 13 AU and r'^; 

G3 - gradient with 40% forsterite at 13 AU and r'^; 

R14 - ring with 60% forsterite between 13 and 14 AU; 

R20 - ring with 40% forsterite between 13 and 20 AU; 

R25 - ring with 30% forsterite between 13 and 25 AU; 

R60 - ring with 10% forsterite between 13 and 60 AU; 

RW30 - only 150 and 200 K forsterite, 60% between 13 and 30 AU; 
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30% require outer radii of 20 (R20) and 25 AU (R25) respec- 
tively. At these radii, midplane temperatures are down to 100 
K (Fig. [5]), and the 69 fim feature is dominated by this compo- 
nent with a peak at 69 //m. In contrast to the gradient models 
Gl and G2, these rings do not contain a 50 K component, since 
these temperatures occur farther out in the disk. Therefore the 
feature shape matches much better than models with a gradient, 
although still some iron is required to shift the feature towards 
longer wavelength. At even lower abundances (R60, 10 % and 
60 AU), the feature strengths at shorter wavelengths cannot be 
reproduced. 

5.2.4. Only warm forsterite 

In order to reproduce the 69 fim feature with iron-free forsterite, 
we need to exclude contributions from 100 K forsterite. 
However, ring models that do so (R14) do not produce enough 
flux at all wavelengths because the emitting region is too small. 
Larger rings include warm forsterite from the disk surface, but 
also colder forsterite from the midplane (R20 and R25). We 
therefore run the same series of ring models, but now with 
forsterite only in regions that are 150 and 200 K, disabling the 
contribution from colder forsterit e. This model is e asily com- 
pared to the optically thin fit of ISturm et all (|20T0.) (Fig. |4]i, 
which also has only these components. We implement this by 
replacing the opacities for 50, 100 and 295 K forsterite with that 
of amorphous silicates. A ring with an outer radius of 30 AU and 
abundance of 60 % provides a good fit to the short wavelength 
features (RW30). The shape of the 69 /vm feature matches the 
observed one well. 

Looking at figure [3] it becomes clear that such a model with 
only warm forsterite has an abundance that varies with height 
above the midplane. Outside of the disk wall (> 15 AU), a ver- 
tical temperature gradient is present from a warm surface to a 
cold midplane. The forsterite in this region must therefore be lo- 
cated only at the disk surface, with an abundance that decreases 
to zero near the disks midplane. The forsterite is present in the 
disk wall and near the disk surface farther out, but not in the mid- 
plane. This scenario ag rees well with the geometry proposed by 
iBouwman et alj (l2003l) . where forsterite is produced in the disk 
gap and transported into the disk wall and outer disk surface. 

5.2.5. Best fit to 69 /zm feature 

The parameter study from the previous section yields a number 
of good fits to the integrated feature strengths, but only 2 models 
fit the shape of the 69 fim feature well: model R20, with 40% 
forsterite between 13 and 20 AU and model RW30, with 60% 
forsterite between 13 and 30 AU, but only in the 'warm' regions 
with a temperature of 150 and 200 K. Fig. |6] shows the mod- 
eled feature shape, as well as the contribution of the individual 
temperature components. 

Model R20 peaks at 69.10 fim, and the model has been 
shifted by 0. 1 ;um to obtain a good fit. We interpret this shift as an 
increase in the iron fraction from to 0.3 %, following the inter- 
polation for low iron content bv lSturm et al.l(l2010i) . which gives 
a peak shift of 3.13 %/iJ.m. The peak shows dominant contribu- 
tions from both 200, 150 and 100 K, and is therefore s lightly 
broader than the optically thin fit from ISturm et al. l (2010). but 
well within the error bars. The higher temperature components 
dominate the feature by integrated emission, but the 100 K com- 
ponent is narrower and is equally important for the peak position. 



Model RW30 peaks at 69.15 fim, and provides a good by- 
eye fit without a wavelength shift. Fit quality increases with an 
additional shift of 0.05 fim, corresponding to an iron fraction of 
less than 0.2 %. The relative contribution of the 200 and 150 
K components are 57 and 42% in contrast with 62 and 27 % 
for the optically thin model (B. Sturm, private communication). 
This also explains the small wavelength shift required, as the 150 
K component has its peak at 69.10 fim. 

The feature of model RW30 is also a little bit too broad, es- 
pecially at the blue side of the feature. The reason for this broad 
feature is the high forsterite abundance of 60%, at which the fea- 
ture itself starts to become optically thick. A lower abundance 
would require a larger region with a temperature of 200 K, which 
is difficult to achieve in the current model. 



5.3. Location of forsterite emission 

To visualize where the different forsterite features originate in 
the disk, we calculate the contribution to the features at every 
location in the disk. Although the temperature and optical depth 
also give some information about where the features originate 
(Fig. [3]), it does not take into account that the density decreases 
with disk height and radius, as well as the increase in opacity in 
the wavelength of the feature and inclination effects. 

To calculate these contributions, the ray tracer stores the in- 
dividual contributions to the total spectrum in every grid cell. 
We then fit a 3rd order polynomial to the continuum, subtract it, 
and calculate the integrated feature strength. The contributions 
of the individual temperature components to the 69 fim feature 
are calculated in a similar fashion. 

Because we view the disk at a 42° angle, the upper half of 
the disk (z/r>0) is in our direct line of sight, while the bottom 
half (z/r<0) is mostly obscured from view. This results in an 
asymmetric profile around the midplane (z/r-Q) in Figures Qand 
|8] while the forsterite abundance itself is symmetric around the 
midplane in all our models. 

5.3.1. Sliort wavelengtli features 

The contributions to the integrated flux of the features at 11.3, 
23.5, 33.5 and 69 fim are plotted in figure Q for best-fit mod- 
els R20 and RW30. All features have a strong contribution from 
the disk wall, approximately half of the total integrated fea- 
ture strength at all wavelengths. The remaining emission orig- 
inates at larger radii at or near the disk surface. Features at 
longer wavelength can be produced by colder material, and have 
stronger contributions from larger radii. Due to the lower contin- 
uum opacity at longer wavelengths, features also originate from 
deeper within the disk, especially at larger radii where surface 
densities become lower. 

The shortest wavelength feature at 11.3 fim originates in the 
warmest regions of the outer disk, mostly in the disk wall. Its 
surface component extends up to only a few AU outwards of 
the disk wall, and is located at the disk optical surface for both 
models. At 24 and 33 fim, the emission from the disk surface 
extends farther out, and its contribution to the total features in- 
creases up to more than 50%. As the optical depth of the con- 
tinuum decreases at these wavelengths, features start originating 
from deeper in the disk. Both models look similar at 24 fim but 
start to deviate from each other at 33 fim, as here the cold 100 K 
forsterite component is not present, truncating most of the emis- 
sion after 17-18 AU. 
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Fig. 6. Modelled 69 ;um feature shape (solid lines) compare to the PACS observations (grey error bars) for models R20 with a 0. 1 
^m shift towards lo nger wavelengths (left) and model RW30 with a 0.05 /vm shift (right). Overplotted are the optically thin fit from 
ISturm etaP (j2010)(grey solid line), and contributions from the individual temperature components. The bottom panels show the 
difference between the PACS obseiTations and the modelled spectra. 



At 69 fim, the optical depth of the continuum is so low that 
features originate from regions all the way down to the mid- 
plane, even in regions just behind the disk wall (14 AU). The 
total contribution from outside the disk wall is more than 50% 
(see also Fig. |9l), and peaks towards the midplane with only a 
very small contribution from the disk optical surface. The model 
with only warm forsterite doesn't extend as far out because the 
cold forsterite component in the midplane is not present. 

The detailed analysis of the feature location confirms the 
simple analysis from the beginning of this section: towards 
longer wavelengths, the optical depth decreases and features 
originate from regions deeper in the disk which are hidden from 
view at shorter wavelengths. Compared to figure |3] regions of 
emission don't reach down as far as the vertical r=l surface. 
The main reason is that disk's inclination is (by definition) not 
taken into account in calculating the vertical t-1 surface. When 
viewing the disk at an inclination angle i, the t= 1 surface of an 
inclined disk is raised by approximately a factor ~l/cos(i) with 
respect to a face-on disk. 

5.3.2. The 69 ^im feature 

The location of emission of the different temperature compo- 
nents of the 69 fim feature are displayed in Fig.|8]for model R20 
and RW30. These are calculated similarly to the shorter wave- 
length features - continuum subtracting the spectrum in every 
grid cell - but for each of the five temperature components sepa- 
rately. For model R20, the cumulative radial contribution to the 
69 yum feature is plotted in figure |9] 

For model R20, the 295 K component originates from the 
disk wall (between 13 and 14 AU), except from near the mid- 
plane. This is because the temperature in the midplane is slightly 
lower than at the surface, an effect of the shadow of the inner 
disk (see Fig. [T]). The 200 K component originates only for a 
small part in the disk wall. Most of the emission comes from the 
upper half of the disk between 13 and 15 AU, and peaks below 



the disk surface. At 150 K, a fraction of the emission comes from 
the shadowed part of the disk wall, and the bulk from near the 
midplane below the surface out to 17 AU. The 100 K component 
does not originate in the disk wall, but from the midplane out- 
wards of 14 AU. The coldest, 50 K component originates from 
even larger radii(> 18At/) in the midplane. 

For Model RW30, the 200 K and 150 K components origi- 
nate in more or less the same regions as for model R20, while 
the 50 K, 100 K and 295 K components are not present. 

6. Discussion 

6. 1 . Forsterite abundance in comparison witli previous work 

The silicate mineralogy of HD 100546 has been addres sed in 
a num ber of studies (e.g. Bouwman et al. 2003, Juha sz et al.l 
1201 Oh . with the underlying assumption that emission features 
arise in an optically thin surface layer in the disk. Assuming 
a homogeneous composition for the emitting region, the mass 
fractio n of forsterite in small du st grains is found to be around 
5-6% (Ivan Boekel et alJi2005 . Juhasz et alj,20 10). Because this 
number can be derived from the mid-infrared spectrum directly 
by measuring crystalline and amorphous feature strengths, we 
will refer to this as the apparent crystallinity. It is clear that 
this number is not neccessarily representative for the mass ra- 
tio in the entire disk, but an extrapolation is of course tempt- 
ing. Taking a separate dust c omposition for the inner and outer 
disk, iBouwman et al. I dlool find a forsterite abundance in the 
depleted inner disk of 2%, compared to 19% in the outer disk. 
Since the latter dominates the total mass, the total mass fraction 
is also 19%, a factor of 4 higher than the apparent crystallinity. 

Our analysis shows that the crystalline forsterite is concen- 
trated towards the disk wall, with locally very high abundances 
(40-60%). Although uncertainties in the grain properties and 
model assumptions can lead to small errors in the derived abun- 
dance, the abundances we find are significantly higher than those 
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Fig. 7. Location of emission of the forsterite features at 11, 
24, 33 and 69 jjm (top to bottom) for models R20 (left) and 
RW30 (right). Overplotted is the radial t=1 surface in the op- 
tical. Contours are given as a fraction of the maximum emission, 
and are between 1 and 100 %. The emission has been integrated 
over the azimuthal direction. The emission locations are asym- 
metric around midplane, because the disk is seen from above, 
and therefore we do not see the bottom half of the disk. Contours 
are drawn as fraction of the peak intensity at 100%, 10% and 1 %, 
and color coded with dark grey between 10% and 100% and light 
grey between 1% and 10% 

Forsterite mass 



Model 


Forsterite mass [Me] 


Fraction [%] 


R20 


0.27 


0.8 


RW30 


0.17 


0.5 


G2 


0.33 


1.0 


5% 


1.7 


5.0 



Table 4. Forsterite mass of our best fit disk models, plus a gra- 
dient and constant abundance model for comparison. 



quoted above. These high abundances are confined to a very nar- 
row region (from 13 AU up to 20-30 AU). Outside of this region - 
where the bulk of the mass is located - the forsterite abundance is 
consistent with zero. Therefore, the overall forsterite mass frac- 
tion - the mass in crystalline forsterite grains compared to the 
total disk mass in small dust grains - is around 0.5-0.8% (See 
Tabled. This is a factor of 10 lower than quoted above. 

The reason why such a small amount of forsterite shows up 
so prominently in the SED is a typical case of window-dressing: 
The wall reprocesses approximately a third of the total stellar 
energy captured by the disk, but represents only a few percent 
of the total disk mass. The wall therefore acts as a display case, 
showing a strong apparent crystallinity from only a small amount 
of forsterite with a high local abundance. This may also explain 
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Fig. 8. Location of emission of the forsterite 69 ;um features, split 
into its different temperature components for models R20 (left) 
and RW30 (right). Overplotted is the radial t-1 surface in the 
optical. Contours are given as a fraction of the maximum emis- 
sion of the total feature, between 1 and 10 % (light grey) and 10 
and 100 % (dark grey). The emission has been integrated over 
the azimuthal direction. 
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the lack of correlations found in larger samples between appar- 
ent crystallinity and other disk parameters. 



6.2. Forsterite abundance in comparison witli Solar System 
comets 

It has been noted before that the mid-infrared spectrum of 
comet Hale-Bopp shows a strik ing resemblance to that of HD 
100546 (e.g. Mal fait et al.lll998h. The derived cry stalUnities are 
in the range of 25 % to 50 % (iBrucato et anil999l; IWooden et all 
11999: Harker et all |2002|) . It should be noted that the degree 
of crystallinity depends on the assumed dust model, see e.g., 
iHanner & Zolenskvl (l2010h . With the inclusion of larger grains 
(>yu m), crystalUnities tend to be lower, aroun d 7% for Hale-Bopp 
(e.g. lMin et a l. 2005b: Bou wman et aDl2003l) . Because we mea- 
sure the crystallinity of small dust grains in HD 100546, we can 
compare the local abundance in the disk wall of 40-60% with the 
crystallinity of small grains in comets. 

Also other longer period comets with high crystallinities 
have been observe d, for example COM ET C/2001 Q4 has a crys- 
tallinity of -30% (IWooden et al.ll2004 )). The surface of shorter 
period comets might be subject to amorphization by solar irradi- 
ation, but ejecta from the inside comet 9P/Tempel released dur- 
ing the Deep Impact experiment sho w that these comets can also 
have a high crystallinity of 20-40% (iHarker et al.ll2007l) . 

Such high crystallinities at these large radii are not expected 
from radial mixing models (Gail 2004), and not observed in typi- 
cal protoplanetary disks. Because the proposed formation region 
of comets is beyond the snowline, it coincides with the location 
of increased crystallinity in the disk wall of HD 100546, making 
it an ideal place for comet formation. 



6.3. Crystallinity-gap correlation 

Both the high abundance and spatial distribution indicate a cor- 
relation between the disk geometry and the strong crystalline 
silicate features in the disk of HD 100546. If the disk gap were 
filled with amorphous material - so that an illuminated wall is 
no longer present - the spectral features of forsterite would be 
weaker by roughly a factor of 3. This might indicate that the 
strength of crystalline features - the apparent crystallinity - does 
not reflect the overall crystal abundance, but rather its spatial 
distribution and the presence of a disk gap. 

Another star which s hows a similar geome try and apparent 
crystallinity is RECX 5 (iBouwman et al.ll201(A . This opens up 
the possibility that a disk gap-crystallinity correlation also exists 
for other protoplanetary disks. However, not all d isks with gaps 
show crystalline dust features (iBrown et al.ll2007 T). A reason for 
this might be that the crystalline features only show up if the gap 
is located in the right temperature range. A good census of crys- 
tallinity and gaps in the right radius/temperature range is needed 
to establish the presence of such a correlation. The strength of 
crystalline emission features ca n be well deter mined from mid- 
infrared spectroscopy (e.g. Juh asz et al.ll2010l) . But the charac- 
terization of disk gaps from the SED is less tri vial, and requires 
spatially resolved observations (IVamiere et al.l l2006). 

Aside from the apparent crystallinity of HD 10 0546 - which 
is one of the highest among Herbig Ae stars (e.g. lJuhasz et al.l 
l2oToh - also the local abundance of forsterite at tens of AU 
is unusually high (40-60% at 13-20 AU). Other Herbig Ae/Be 
stars show lower forsterite abundances in the outer (2-20 AU) 
disk, around 10-15% d erived with an optically thin approach 
dvan Boekel et all l2004i) . with the exception of HD 142527, 



which also has a prominent disk gap (IVerhoeff et alJl20Tll) . If 
all Herbig Ae/Be stars would have the same crystallinity as HD 
100546 in the 13-20 AU range - and an equal or higher abun- 
dance farther in - they would produce much stronger spectral 
features than observed. Only if the abundance would increase 
outwards would it be possible to produce the observed crys- 
tallinities in the 0-5% range with the same high local abun- 
dance at 13 AU. However, this is in sharp contrast with o bser- 
vations (jv an Boek el et a n i200l and radial mixing models 
l200lllGaill200 4). that show a crystalUnity that is constant or de- 
creases outwards. 

Whether the spatial distribution of forsterite around HD 
100546 is truly unique among Herbig Ae/Be stars requires an 
analysis of disks with lower crystallinity as well. The Herschel 
Open Time Key Program DIGIT will observe the 69^m feature 
in a number of Herbig Ae/Be stars, allowing for a characteri- 
zation of the spatial distribution of forsterite in these disks in 
a similar fashion. If there is a mechanism that produces a high 
crystallinity at 13-20 AU independent of a disk gap, it will be 
identified in this sample. 

6.4. Forsterite origin 

How can such a high abundance of forsterite originate in the disk 
wall, and is it related to the formation or presence of the disk 
gap? There are two chemical pathways of producing forsterite: 
by condensation from the gas phase or by annealing from amor- 
phous olivines. Both require high temperatures above the dust 
evaporation temperature (~ 1500 K) or the glass temperature (~ 
1000 K) respectively. These high dust temperatures exist only 
in the inner regions of protoplanetary disks, but not at 13AU. 
We can think of two main scenarios, radial transport or in-situ 
production: 

In the first scenario, forsterite is produced in the inner disk, 
and transported outwards to 13 AU. The current dust mass in the 
inner disk is much lower than the total amount of forsterite in 
the disk wall between 13 and 20 AU. It therefore seems unlikely 
that radial mixing from the current inner disk can produce such 
a high crystallinity in the disk wall. As discussed before, it is un- 
likely that radial mixing has produced a locally high abundance 
before the opening of the disk gap, as it would have shown up in 
the mid-infrared spectra. This means that radial transport must 
have taken place during the formation of the disk gap and deple- 
tion of the inner disk. Radiation pressure on the crystalline grains 
can transport them across the gap and blow them into the disk 
wall, depleting the inner disk. Also a different mechanism, pho- 
tophoresis, can tran sport grain s outwards efficiently when the 
inner disk depletes (iPetit et al.l 12006). The depleted dust mass 
from the inner disk (200 times the current inner disk mass, sec- 
tion l4.3b is large enough to explain the crystalline silicates in the 
disk wall, but it requires extremely efficient radial transport and 
almost a 100% crystallinity in the inner 4 AU before the opening 
of the gap. 

In the second scenario, forsterite is produced locally in the 
disk wall. Because the equilibrium temperature of the dust in the 
wall (200 K) is far below the dust evaporation or glass temper- 
ature, another mechanism must then be responsible for heating 
the amorphous dust grains so they can crystallize: 



- Flash heating of grains in a stellar outburst ([Abraham et al.l 
2009). However, HD 100546 shows currently no signs of 
such outbursts. Such crystals are observed to disappear 
quickly after formation, which is not the case for HD 
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100546, where crystals persist over more than a decade of 
observations. 

- Shocks in the disk midplan e can locally heat the grains 
above the glass temperature (iHarker & DeschI [20021) . From 
scattered light images, HD 100 546 is known to possess 
two spiral arms in the outer disk (iGradv et al.ll200ll) . which 
might be related to shocks in the disk midplane that could 
create forsterite (Harker & Desch 2002). If a companion is 
responsible for clearing the disk gap, it likely induces a 
similar spiral wave at the location of the disk wall (e.g. 
iNelson & Papaloizoij2003l) . producing forsterite locally. 

- Amorphous silicates can be crystallized in parent bodies 
large enough to retain the heat from decay of radioactive 
elements or accretion (e.g. iHuss et aTl 1200 Ih . Micron-size 
forsterite grains can be created in a collisional cascade of 
such objects, tha t can take place in the disk gap or wall (see 
iBouwman et al . 2003 for a discussion on this topic.) The 
amount of forsterite found in the disk wall of HD 100546 re- 
quires the destruction of an Earth-Uke planet or several hun- 
dred Pluto-sized objects. 

All plausible scenarios for forsterite formation (shock waves, 
parent body processing and radial transport during gap forma- 
tion) are connected to the presence or formation of the disk 
gap. Therefore we conclude that the anomalously high forsterite 
abundance in the disk wall of HD 100546 is - one way or the 
other - connected to the presence of the gap. 

6.5. Forsterite in the inner disk 

The mass fraction of forsterite in the inner disk is hard to con- 
strain from our models. The contribution to mid-infrared fea- 
tures in the model is low, and shows a significant contribution 
only at 1 1 jjm. Regions near the inner rim that are above a tem- 
perature of 1000 K should be fully crystalline, but even if the 
silicates in the entire inner disk are fully crystalline, the contri- 
bution is only a few Jy fim at 1 1 yum. The contribution to longer 
wavelengths is negligible because the inner disk is too warm, and 
optical depth effects do not play a role b ecause it is also (verti- 
cally) optically thin (Benist v et al.ll20 10'). Therefore we cannot 
put any limits on the crystallinity of the inner disk with the cur- 
rent model and observations. 

6.6. Tliermal contact 

We assume that the different dust species in our model (amor- 
phous silicat es, forsterite , carb on) are in thermal contact. 
However, Bouwm an et a D (2003) found with an optically thin 
approximation that the bulk of the forsterite is at a temperature 
of ~ 70 K while the other dust species are much warmer, leading 
them to the conclusion that forsterite cannot be in thermal con- 
tact with the other dust species. Because the 2D radiative transfer 
does not need this cold component due to optical depth effects 
(see section 15.1b . the forsterite has the same temperature as the 
other dust species in the disk wall, which means they can be in 
thermal contact. 

If the forsterite in the outer disk is forced to be out of thermal 
contact with the other dust species, it does not reach tempera- 
tures high enough to explain the shortest wavelength features of 
forsterite. The reason for this is the very low opacity of forsterite 
at optical wavelengths compared to the other dust species, es- 
pecially carbon. Therefore the forsterite is not heated directly 
by the central star, but only by mid-infrared radiation from the 
disk surface. It will therefore have a temperature similar to dust 



species in the midplane (50-100 K), which is too cold to explain 
any short-wavelength features. Even if the forsterite is moved in 
to the inner disk, where it would normally be hotter, it falls short 
of explaining the strength of the short wavelength features (see 
also section 1631) . 

We can ther efore exclude th a t the dust species are not in 
thermal contact. Bouwman et al.l (1200 3) argue that the dust in 
HD 100546 is more primitive than in comet Hale-Bopp because 
it is out of thermal contact, and that we see the creation of such 
comets rather than their destruction. Since we exclude that the 
dust species are not in thermal contact, this implies that the 
dust in HD 100546 is more similar to comet Hale-Bopp than 
thought before, and not necessarily more primitive. This opens 
up again the possibility that we are witnessing the destruction of 
forsterite-rich comets like Hale-Bopp in HD 100546, rather than 
their formation. 

6. 7. Iron content 

Finally, we use the shift in wavelength required to fit the shape of 
the 69 yum feature to estimate the iron content of the silicates. For 
large changes in the iron fr action, the central wavelength shifts 
linearly with iron content (IKoike et al.ll2003|) . The wavelength 
shifts between an d 10% have not bee n measured, but relies on 
the interpolation bv lSturm et al.l(l2010l) . which assumes the shift 
in this regime is also linear 

We also assume that the shape of the 69 fim feature does 
not change with a small change in iron content. Although this 
has not been measured in the lab for forsterite with low iron 
fractions, small changes in iron fraction at higher iron content 
have been measured ( Koike et al. 20 03). For example, between 
9.3 and 8.6%, the shape does not change significantly. The fea- 
ture broadens a little bit between and 8% iron content, but this 
change is much smaller than the broadening due to a temperature 
change (see Fig. lC.ll ). 

Therefore we expect the shape of the feature not to change 
significantly for the small changes in iron content (<0.3%) found 
in this paper. Laboratory measurements of forsterite with low 
iron content are necessary to confirm this. 

7. Conclusion 

We have studied the spatial distribution of crystalline silicate 
dust (forsterite) in the gapped disk of Herbig Be star HD 100546. 
We have used a 2D radiative transfer code to model the dust ge- 
ometry, and study the effects of optical depth on the forsterite 
feature strengths in the mid and far-infrared, as well as the 
shape of the temperature-dependent 69 fim feature observed with 
Herschel PACS. Our conclusions are: 

1. 2D radiative transfer effects need to be taken into account 
when interpreting spectral features over a broad wavelength 
range. In the case of HD 100546, the bulk of the warm (200- 
150 K) forsterite that contributes to the 69 yum feature is hid- 
den from view at shorter wavelengths (10-40 fim). An addi- 
tional cold forsterite component (~50-70 K) is therefore not 
necessary to explain the strength of the 69 fim feature. 

2. The forsterite of HD 100546 is located in the disk wall at 
13 AU, with extremely high local abundance. Based on mid- 
infrared features strengths and the shape of the 69 jum fea- 
ture, we arrive at two best-fit models: 

a) a model with 40% forsterite between 13 and 20 AU. 

b) a model with 60% forsterite between 13 and 30 AU but 
restricted to regions where the temperatures are around 150 
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and 200 K. This excludes the cold midplane outside of the 
disk waU. 

3. All forsterite features show a strong contribution from the 
disk wall, but also probe a region farther out. The 69 fim fea- 
ture from this region originates from below the disk surface. 
The shorter wavelength features from this region originate at 
the disk surface. 

4. The mass fraction of forsterite compared to the total disk 
mass is low - around 0.5-0.8% - a factor of 10 lower than the 
apparent crystallinity derived from and optically thin analy- 
sis of relative feature strengths in the mid-infrared. The wall 
at 13 AU acts as a display case for showing a high local 
forsterite abundance. 

5. The crystalline forsterite in the disk of HD 100546 is in ther- 
mal contact with the other dust species. 

6. We cannot rule out forsterite that is completely iron free. Our 
best fits indicate an iron content of <0.3 % of the crystalline 
olivines. 

7. Due to a degeneracy in the modelling between iron content 
and dust temperature, it is not straightforward to use the 69 
jum feature as an independent temperature indicator. 

8. A disk model in hydrostatic equilibrium (with Tgas=Tdust) 
can self-consistently explain the SED and mid-infrared 
imaging. Our best fit has a surface density power law oc 
and an inner disk depletion factor of 200. Additionally, con- 
tinuum opacity sources - such as carbon - must have a low 
abundance to fit the low flux at 8 /im. 

9. The apparent crystalhnity of dust as measured from the ra- 
tio of crystalline and amorphous components of the mid- 
infrared features is not necessarily a good measure of the 
real crystallinity, the mass fraction of small dust grains in 
the disk that is crystalline. Both disk geometry and spatial 
distribution of the crystals can significantly distort the pic- 
ture, and a one-to-one correlation between crystal feature 
strength and disk state is not to be expected. However, in our 
case, the detailed analysis shows a highly localized region 
with high abundance, which may very well inform us about 
special processing in this region, possibly linked to planet 
formation. 
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Fig. A.l. VISIR Q band image at 18.7 //m. North is up and East 
is left, with a 3x3 "field of view. The scale is logarithmic in 
intensity. 
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Appendix A: VISIR imaging 

AA. VISIR observations 

During the night of March 26 2005, Q-band imaging was per- 
formed usi ng the VLT Image r and Spectrometer for the mid- 
IR (VISIR. iLagage et alj|2004l) . For photometric calibration and 
PSF determination, the standard stars HD 50310 and HD 150798 
were observed before and after the science observation respec- 
tively. Standard chopping and nodding was employed to remove 
the atmospheric background emission. The imaging was per- 
formed with the Q2 filter {A^ = 18.72 yum) in the small field mode 
(pixel field of view - 0.075"). The observing conditions were 
fair, airmass ~1.4, and the optical seeing ~0.8". The achieved 
sensitivity was 80 before and ~100mJy/10crlh after the science 
observation. The full width at half maximum (FWHM) of the 
point spread function (PSF) was -0.5". The data reduction was 
performed with a dedicated pipeline, which corrects for vari- 
ous in strumental signatures (see Paiitin et al. 2008, 2009; Pantin. 
120101) . The reduced image is presented in figure I A. l1 

A.2. Image analysis 

The VISIR image (fig lA.ll) shows a resolved disk that is almost 
completely spherical, consi stent with previous imaging at this 
wavelength (i Liu et al]|2003h . Deviations from spherical symme- 
try are down to the 10 % level, which means that inclination 
must be small. Due to the relatively large size of the PSF at this 
wavelengths, it is difficult to get a direct constraint on the in- 
clination out of the image. Comparison with a set of inclined, 
PSF-convolved model images show that it is consistent with an 
inclination less than ~40 degrees and a position angle along the 
South-East to North- West. 

Due to the increased sensitivity with respect to previous ob- 
servations, the image shows resolved emission up to ~1.4", cor- 
responding to a physical radius of 145 AU at 103 pc. To char- 
acterize this emission, we construct a surface brightness pro- 
file. Because the disk is not seen face-on, we construct the sur- 
face brightness profile from the image by averaging the surface 
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Fig. A.2. Radial emission profiles in the Q band at 18.7 /im, cal- 
culated assuming elliptic annuli (see text for details). Displayed 
are: VISIR image (diamonds) and reference PSF (grey line). The 
errors on the PSF are displayed as horizontal error bars, those on 
the noise as vertical ones. Overplotted are PSF-convolved disk 
models with a varying surface density power law of p=0.0, 1.0 
and 2.0 (dotted, solid and dashed respectively). NOTE: radius 
refers to the semi-major axis of the ellipse, which reflects the 
physical radius where the radiation is emitted. 



brightness over an elliptical annulus - rather than a circular one. 
The shape of the annuli correspond to a disk inclination of 42° 
and a position angle of 145°. The semi-major axis of this ellipse 
corresponds to the distance from the star at which the radiation 
is emitted (radius). Annuli are 0.075" wide, the size of one pixel 
on the VISIR chip. They are thus much smaller than the width 
of the VISIR PSF 

Comparing the measured radial intensity profiles to that of 
the VISIR PSiQ (Fig lA!2l i, we see that the disk is clearly re- 
solved. The profile itself is not easily characterized by simple 
power-laws as is often the case for scattered light images, due 
to the different emission mechanism (thermal versus scattered 
light) as well the relatively large size of the PSF. 

Instead, we recognize two components. The profile up to 
0.5" is well described by a Gaussian with a FWHM of 0.73". 
After quadratic subtraction of the PSF, we derive a spatial scale 
of 0.33" for this central component, corresponding to a physical 
scale of FWHM - 34 AU at a distance of 103 pc. This num- 
ber agrees well with the value of 34 + 2 found by ILiu et^ 
(2003) with direct imaging at the same wavelength. This com- 
ponent therefore originates in a region just behind the disk wall, 
which has a spatial scale of ~ 26 AU. 

The second component - outside of 0.5"- comes from farther 
out in the disk, and has not been analysed before at this wave- 
length. Extended emission is seen up to 1.4", corresponding to 
a physical radius of 145 AU at 103 pc. This emission probes the 
thermal continuum emission from the outer disk surface. The 
bump around 1.0" corresponds to the first diffraction ring in the 
VISIR PSF, and does not represent a real structure in the disk. 
Because the radial profile traces the disk surface over almost an 
order of magnitude in radius, we can use it to put limits on the 
surface density of small grains, which we will do in the next 
section. 

^ To be consistent, the plotted surface brightness profile of the PSF 
is also calculated using an elliptical annulus even though it is spherical, 
and its FWHM in the plot is larger than the 0.5"previously quoted. 
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Fig. B.l. Example of continuum subtraction, for model R20. The 
top panel shows the model spectra (black line) and calculated 
continuum (grey line). The bottom panel shows the continuum 
subtracted spectra, with the same color-coding. 
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Fig. B.2. Continuum subtraction for the Spitzer spectrum. The 
top panel shows the model spectra (black line) and calculated 
continuum (grey line). The bottom panel shows the continuum 
subtracted spectra, with the same color-coding. 

A.3. limits on the surface density distribution of small grains 

To test if the assumed power law index of p=l for the surface 
density profile matches our VISIR image, we compare model 
images of varying SDP with the observed image, keeping the 
dust mass and outer radius fixed, as well as the inner disk surface 
density. Because the pixel size is much smaller than the VISIR 
PSF, we convolve it with our model images before constructing 
the radial intensity profiles, and take into account the errors on 
the width of PSF (FWHM = 0.509±0.008") in the fitting proce- 
dure. 

Because we solve for the hydrostatic disk structure, the shape 
of the radial profile depends only on the surface density distri- 
bution (Figure lA!2l i. The model is best fit for a surface density 
power law with index p 
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Appendix B: continuum subtraction 

To measure feature strengths and compare continuum subtracted 
spectra, we have to fit and subtract the continuum from both our 



model spectra and data. An example for model R20 is shown in 
figure IB. II for the Spitzer spectrum in figure IB. 21 We are inter- 
ested in the narrow emission features of forsterite (a crystalline 
silicate) which are superimposed on much broader features from 
amorphous silicates, that occupy the same wavelength region - 
except for the 69 yum feature. The continuum - that includes these 
broad features - is therefore not easily described by a low-order 
polynomial. 

The continuum is better fit using a spline function, which is 
essentially a piecewise polynomial between (arbitrary) contin- 
uum points - or knots. At every knot, the polynomials are con- 
nected such that they have a continuus derivative. The best result 
is obtained when knots are placed just outside the features, and 
the knots ('continuum points') between features are spaced by a 
distance roughly equal to the feature width. 

The fitted spline follows the shape of the continuum out- 
side the features very well. Inside the features, the continuum 
is smooth, but has more structure than to be expected from the 
opacities of amorphous silicates only. Especially the bumps at 
12, 27 and 32 fj.m are artifacts from the continuum subtraction, 
and not from the underlying amorphous silicate features. They 
arise because the opacity of forsterite between features is not 
always zero, especially between the big features and the little 
'wings' at 11, 24/28 and 33 //m. These wings are much weaker or 
not seen in the Spitzer observations (Fig.|4]i, and do appear much 
weaker in other forsterite opacities (Servoiri & Piripu 1973). 

When comparing the model to the data, we therefore want 
to exclude these wings, and choose our knots between feature 
and wings. Although the resulting continuum is not completely 
smooth - which systematically underestimates the integrated flux 
- we apply the same continuum subtraction to the observations. 
This cancels out the systematic error, allowing for a comparison 
of the shape of the feature, without a systematic error from the 
wings - which is not in the data. For the 33 /urn feature, there is 
the additional complication that the Spitzer spectrum runs only 
up to 36.9 jum, and it is not possible to choose a continuum 
point outside the little wing. We therefore force the tangent at 
this point to zero to get a reasonable continuum, although this 
makes the 33 jum flux appear a little stronger than it really is. 



Appendix C: Implementation of 
temperature-dependent opacities 

When we implement temperature-dependent opacities into the 
radiative transfer code, we have to take into account that the 
opacities that set the temperature depend on that same tempera- 
tur^. We therefore recalculate the opacities at every iteration of 
the hydrostatic vertical structure loop, to keep the solution self- 
consistent. 

Optical constants for forsterite have been measu red at tem- 
peratures of T,= 50, 100, 150, 200 and 295 K jSuto et alj 
2006) for iron free forsterite (Mg2Si04), but no measurements 
exist for crystalline olivines with small iron contaminations 
(FejMg2-ASi04). For temperatures between 50 and 295 K, we 
calculate the opacity from the optical constants for every tem- 
perature component separately, and combine them using a linear 
interpolation. The opacity at temperature T is calculated as: 



K(T) 



Ti 



- i+\ 



Ki+l 



(C.l) 



' Although in the case of forsterite, the changes in the opacities are 
too small to have a major influence on the disk structure 
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Temperature dependent forsterite opacity 
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Fig. C.l. Temperature dependence of the forsterite opacity 
around 69 fim, calculated for the dust grains used in this paper 

where /c,- is the opacity at temperature T,, and T is between T,- and 
T,+i . Below T= 50 K, we use the opacities of 50 K forsterite, and 
we use the 295 K opacity above a temperature of 295 K. This 
approach does not take into account the shift in wavelength that 
might occur outside the 50-295 K range. However the pe ak shift 
becomes smaller at lower temperatures feuto e t aTll2006l) and is 
not so important below 50 K. Temperatures above 295 K are not 
present in the SED due to the disk gap, except in the inner disk 
which is does not contribute significantly at 69 fim because it is 
much hotter 



